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Context. The origin of subluminous Type la supernovae (SNe la) has long eluded any explanation, as all Chandrasekhar-mass models 
have severe problems reproducing them. Recently, it has been proposed that violent mergers of two white dwarfs of 0.9M o could lead 
to subluminous SNe la events that resemble 1991bg-like SNe la. 

Aims. Here we investigate whether this scenario still works for mergers of two white dwarfs with a mass ratio smaller than one. We 
aim to determine the range of mass ratios for which a detonation still forms during the merger, as only those events will lead to a 
SN la. This range is an important ingredient for population synthesis and one decisive point to judge the viability of the scenario. In 
addition, we perform a resolution study of one of the models. Finally we discuss the connection between violent white dwarf mergers 
with a primary mass of 0.9M o and 1991bg-like SNe la. 

Methods. The latest version of the smoothed particle hydrodynamics code Gadget3 is used to evolve binary systems with different 
mass ratios until they merge. We analyze the result and look for hot spots in which detonations can form. 

Results. We show that mergers of two white dwarfs with a primary white dwarf mass of « 0.9 M Q and a mass ratio larger than about 
0.8 robustly reach the conditions we require to ignite a detonation and thus produce thermonuclear explosions during the merger 
itself. We also find that while our simulations do not yet completely resolve the hot spots, increasing the resolution leads to conditions 
that are even more likely to ignite detonations. Additionally we compare the abundance structure of the ejecta of the thermonuclear 
explosion of two merged white dwarfs with data inferred from observations of a 1991bg-like SN la (SN 2005M). The abundance 
distributions of intermediate mass and iron group elements in velocity space agree qualitatively, and our model reproduces the lack of 
material at high velocities inferred from the observations. 

Conclusions. The violent merger scenario constitutes a robust possibility for two merging white dwarfs to produce a thermonuclear 
explosion. Mergers with a primary white dwarf mass of a 0.9 M Q are very promising candidates for explaining subluminous SNe la. 
This would imply that subluminous SNe la form a distinct class of objects, that are not produced in the standard single white dwarf 
scenario for SNe la, but instead arise from a different progenitor channel and explosion mechanism. 

Key words, stars: supernovae: general - hydrodynamics - binaries: close 



1. Introduction 

Type la supernovae are among the brightest objects in the 
Universe. Their apparent homogeneity makes them one of the 
most important probes for cosmic expansion. This application 
relies on an empirical luminosity calibration (e.g. Phillips 1 993] >. 
Only recently this calibration has been justified by theoretical 
studies, as models were able to reproduce it in terms of the un- 
derlying physics (|Kasen et al.12009 1 |Mazzali etak]|2001) . The 
standard SNe la ( |Branch et al.|l993 1 used as cosmological dis- 
tance indicators make up about 70% ( Li et al.|2010) of the ob- 
served SNe la. There is, however, a number of subclasses of SNe 
la consisting of very dim (e.g. SN 19 91bg, [Le ibundgut et al. 
1993) , normal, bright (e.g. SN 199 IT, |Hhppenko et al.J|l"992J 
and very bright supernovae (e.g. SN 2007if, Sc alzo et al.|2 010). 
Until today, neither the progenitor scenario nor the explosion 
mechanism of SNe la are completely understood. In particular, 
very dim and very bright objects are difficult to explain. There is 
general consensus only that all SNe la are thermonuclear explo- 
sions of white dwarfs (WDs) in binary systems. The nature of 
the companion star and the details of the explosion mechanism, 
however, are still debated. 



Current progenitor models distinguish between the single de- 
generate and the double degenerate scenario. The former (pro- 
posed by Whelan & Iben 1973] > assumes a main sequence or 
giant star, i.e. a non-degenerate star, to be the companion of 
the exploding white dwarf. It feeds the white dwarf via Roche- 
lobe overflow or by winds until the white dwarf approaches the 
Chandrasekhar mass and explodes. In contrast, the double de- 
generate scenario ( |Iben & Tutukov|[T9 84; Webbink 1984) 



as- 



sumes a binary system of two carbon-oxygen white dwarfs. 
Despite being favored by stellar population synthesis studies 



(e.g. Ruiter et al. 2009 1, this scenario has received little atten- 
tion (see, e.g.|Hillebran di~& Niemeyer|2000| i. |Benz et al.| (l990) 
found that in a binary system of a 1 .2 M and a 0.9 M Q white 
dwarf the less massive white dwarf is disrupted and accreted 
onto the more massive white dwarf. This happens within a few 
orbits after mass transfer between the stars has started. It has also 
been shown (e.g. |Nomoto & rben|1985||Saio & Nomoto||1998| l 
that the most likely result of such a merger is an inward propagat- 
ing deflagration flame that converts the white dwarf into an O-Ne 
white dwarf. This white dwarf then undergoes accretion-induced 



collapse rather than thermonuclear explosion (Nomoto & Kondo 
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|1991| . Recent high-resolution simulations of white dwarf merg- 
ers with careful treatment of the initial conditions ( jD'Souza et al.| 
[20061 |Motl et al.||2007) |Dan et al.|[2008l >, however, show indi- 
cations that the immediate disruption of the less massive white 
dwarf may be avoided. Instead, they find stable mass transfer 
for more than 50 orbits. When the less massive secondary WD 
is disrupted after these orbits, the accretion onto the primary is 



slower and carbon shell ignition may be avoided ( |Yoon et"aL 
[2007). 

In the Chandrasekhar-mass scenario, a white dwarf is 
thought to ignite near the center. At first the flame propagates 
subsonically as a deflagration, but burns only a fraction of the 
star. The energy release expands the star and in a later stage of 
the explosion, a detonation triggers that propagates supersoni- 
cally. The mechanism that leads to the formation of the detona- 
tion is still an open question. One of the most promising scenar- 
ios is a deflagration to detonation transition by turbulence (e.g. 
Kho khlov et"aT1[T997| |R6pke & Niemeyer|[2507l |Ropke et al. 



2007a i. Other proposed mechanisms are th e gravitationall y con- 
fined detonation (e.g. |Jordan et al.||2008[ |Meakin et al.||2009| l 



and pulsating reverse detonation (Bravo & Garcfa-Senz 2006) 
Note that a thermonuclear explosion may also be possible before 
the white dwarf reaches the Chandrasekhar-mass if the carbon- 
oxygen white dwarf accretes helium from a non-degenerate or 
white dwarf companion ( |Fink et al.|2007l[20T0)|Sim et al.|20"T0 
IGuillochon et al.|2010] >. 

A more detailed comparison of different explosion mecha- 
nisms can be found in Hillebrandt & Niemeyer (2000 ). However, 
none of the Chandrasekhar-mass models are able to explain very 
dim SN la with a 56 Ni mass of 0.2 M or less. 



Recently, Pak mor et al.| ( |2010| ) proposed that these SNe la 
are the result of a violent merger of two equally massive white 
dwarfs. Once a binary system of two white dwarfs has formed, 
gravitational wave emission leads to continuous shrinking of its 
orbit. At some point, the two white dwarfs are close enough 
to interact directly. If the mass ratio of the two white dwarfs 
is close to unity, the interaction becomes very violent and the 
slightly less massive, and therefore less compact white dwarf is 
destroyed. Its material plunges onto the remaining white dwarf, 
is compressed and heats up. At the interface between the merg- 
ing WDs material can become hot enough for carbon to ignite. 
As carbon burning starts, it is heated up further. If these hot spots 
occur at high enough densities, they can reach the conditions that 
are necessary to ignite a detonation. 

The detonation then propagates through the merged object 
and burns most of the original carbon-oxygen material to heav- 
ier nuclei. As the detonation flame moves supersonically with 
respect to the fuel and the timescale of the flame to cross the 
whole star is much shorter than dynamical timescales involved, 
the system is burned nearly instantaneously at the moment the 
detonation is ignited. With enough energy released from nu- 
clear burning to overcome the gravitational binding energy, the 
merged object explodes, and finally reaches a state of homolo- 
gous expansion. 

In this work, we explore mergers of white dwarfs with a pri- 
mary mass of 0.9 M Q and mass ratios between 0.8 and one. We 
investigate how the conditions during the merging event change 
with mass ratio and determine for which range of mass ratios the 
merger will produce a thermonuclear explosion. 

The paper is structured as follows. In Section [2j the SPH 
code used in this work is summarized and the modifications are 
presented that are introduced to adapt the code to the problem. 
The initial setup is described in Section[3] Section|4]investigates 
inspiral and merger of one specific binary system and discusses 



the conditions that lead to the ignition of a detonation. Section 
[5]compares mergers with different mass ratios. In Section[6]we 
compare the outcome of a thermonuclear explosion following 
a violent merger with observations. Conclusions are drawn in 
Section [7] 



2. Implementation 

In order to simulate the inspiral of the binary system, the lat- 
est version of the smoothed particle hydrodynamics (SPH) code 
Gadget (Springel 2005) is used with some modifications. These 
are needed as Gadget was originally written for cosmological 
applications. The most important modifications to the code are 
the implementation of a stellar equation of state and a nuclear 
reaction network. For the equation of state (EoS) of the elec- 
tron gas, positrons and radiation we use a tabulated version of 

1999[), in which we interpolate 



the EoS by Timmes & Arnett 



linearly. The ions are treated as a fully ionized ideal gas. 

We implemented a 13 isotope nuclear reaction network. The 
network contains all a-elements up to 56 Ni ( 4 He, 12 C, ie O, 20 Ne, 
24 Mg, 28 Si, 32 S, 36 Ar, 4l) Ca, 44 Ti, 48 Cr, 52 Fe, 56 Ni ). The nuclear 
reaction rates are taken from the latest (2009) release of the 
REACLIB database ( |Rauscher & Thielemannp OOO). The net- 
work is solved in each timestep for all active particles whose 
temperatures exceed 10 6 K. It is solved separately from the hy- 
drodynamical equations employing operator splitting. Typically 
the network requires significantly smaller timesteps than the hy- 
drodynamical timesteps; therefore multiple network timesteps 
have to be performed during one hydrodynamical timestep. The 
change in composition is taken into account for the EoS. The en- 
ergy release (or consumption) due to changes in the composition 
is added to the internal energy of a particle. As a consequence 
of the implementation of the stellar equation of state and the 
nuclear network we solve the energy equation instead of the en- 
tropy equation (as in the original version of Gadget). Details and 
tests of these modifications will be discussed in Pakmor ( 201 1) >. 

The configuration of the Gadget code is very similar to its 
previous use in Pakmor et al. (2008). All particles have the same 
mass. The smoothing length is chosen such that each particle has 
50 neighbor particles. To calculate the gravitational forces only 
the tree method is used. The gravitational softening length of a 
particle is set equal to its smoothing length. 



3. Simulation setup 

In a first step, a one-dimensional model of a carbon-oxygen 
white dwarf in hydrostatic equilibrium is constructed. This is 
done by integrating the equation of hydrostatic equilibrium over 
the mass profile of a white dwarf, starting at the center with a 
given central density and assuming a constant temperature and 
nuclear composition throughout the whole star. The central den- 
sity is chosen to give the desired mass of the white dwarf. 

All white dwarfs are assumed to have a constant temperature 
of T = 5 x 10 5 K and a uniform composition of Xq — Xq = 
0.5 is used. To set up the binaries we create white dwarfs of 
M = 0.90 M Q , M = 0.70 M Q , M = 0.76 M G , M = 0.81 M and 
M = 0.89 M , respectively. 

In a second step, each one-dimensional white dwarf model 
is mapped to a distribution of SPH particles resembling its one- 
dimensional profile. Details of this procedure will be described 
in Pakmor (201 1 ). All particles of all white dwarfs represent the 



same mass. The particle mass is chosen such that the most mas- 
sive white dwarf (M = 0.90 M Q ) consists of 10 6 particles. Using 
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Fig. 1. Snapshots of the evolution of the binary system during the inspiral. The system contains two white dwarfs of 0.90 M Q and 
0.81 Mo, respectively. It was set up with an initial orbital period of 33s, rotating clockwise. Color-coded is the logarithm of the 
density on the left and of the temperature on the right. 



the Gadget code, every model is finally relaxed for 100 s to damp 
out spurious numerical noise introduced by the setup. 

After relaxation, the white dwarf models are used to con- 
struct binaries. Both stars in such a binary are set onto a circular 
orbit around their combined center of mass. All particles of each 
star are given the same initial velocity. The period is chosen to 
provide a marginally stable orbit, i.e. the system is stable for 
more than one complete orbit before tidal interactions lead to 
dynamical interactions between the white dwarfs. 

Note that this setup does not relax the binary system once the 
white dwarfs feel the tidal forces of their companion. Performing 
such an additional relaxation step may lead to more stable initial 
conditions and therefore a less violent merger, but to investigate 
this question would be computationally very demanding. 



4. Inspiral and merger of a 0.90M o + 0.81M o white 
dwarf binary 

The binary evolution in the inspiral and merger phase is illus- 
trated in Fig. [T] for the example of a merger between a 0.90 M Q 
and a 0.81 M white dwarf. Shown are snapshots of density and 
temperature in slices through the orbital plane (z = 0). The bi- 
nary system rotates clockwise. It is stable for about two orbits, 
but both stars are slightly distorted by tidal forces. In the begin- 
ning, only the less massive white dwarf is distorted, while the 
more massive star stays relatively unaffected for much longer. 

Mass transfer sets in a few seconds after the start of the simu- 
lation. After about two orbits 60 s) the secondary white dwarf 
starts to break up and its material is accreted rapidly onto the pri- 
mary white dwarf. In this process, material at the surface of the 
primary white dwarf is strongly heated up, first by compression 
and later also by nuclear reactions to temperatures well above 
10 9 K. This process starts with compression that subsequently 
leads to the ignition of carbon burning. The temperature increase 
beyond 1.5 x 10 9 K is due to the energy release in this nuclear 
reaction. 



The evolution of the binary is qualitatively similar to previ- 
ous studies of merging double degenerate systems (e.g.[i 



et al. 2004 Yoon et al. 2007 Loren-Aguilar et al. 2009|>. In these 



studies, the less massive WD is disrupted in the merger and its 
material builds up a hot envelope around the remaining white 
dwarf. Loren-Aguilar et al. ( 2009 ) also find that when material of 
the secondary first hits the primary, the temperature at the inter- 
face becomes high enough to start carbon fusion. In their simula- 
tions, this leads to an rapid increase of the temperature there due 
to nuclear energy release. As the temperature increases, the ma- 
terial becomes non-degenerate, expands and cools down, thereby 
quenching the nuclear reactions. In the end, only a very small 
amount of material is burned. 

There is, however, an important difference in the interpreta- 
tion. Loren-Aguilar et al. ( 2009| l do not consider the possibility 
of a detonation forming at the hot-spot. In this case the detona- 
tion will be able to burn the whole white dwarf. This is a direct 
consequence of the fact that the detonation propagates superson- 
ically and therefore burns all the material before it can expand. 

Hence, the crucial question that decides on the fate of a vi- 
olent merger event is whether the merging process leads to the 
formation of a detonation or not. Unfortunately, the conditions 
and the process leading to the formation of a detonation are still 
not completely understood, but a topic of active research (e.g. 
|Seitenzahl et al.|2009||Ropke et al.|2007b| . 

4.1. Detonation ignition mechanisms 

There are basically two ways a detonation can trigger: Firstly, 
a strong enough shock may lead to direct formation of a deto- 
nation. In this case, an already existing shock wave compresses 
and heats the material causing nuclear burning behind the shock. 
If the nuclear energy release behind the shock becomes large 
enough, it may be sufficient to create a self-sustained detona- 
tion (e.g. Body 1997 1. As the merger is very violent, we nat- 
urally expect strong shocks to occur in the interaction region. 
Unfortunately, however, we are far from being able to resolve 
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Fig. 2. Density vs. temperature scatter plot for all particles of the merger of a 0.90 Mo and a 0.81 M Q white dwarf. 



these shocks on scales small enough to observe the formation 
of a detonation. As all SPH particles in our simulation have the 
same mass, the mass resolution is constant. Therefore regions 
with higher density are resolved better than regions with low 
density. Given that the interaction between the two white dwarfs 
occurs at comparatively low densities, their spatial resolution in 
the interaction region is intrinsically worse than at the center of 
the white dwarfs. In fact, we only barely resolve shocks in the 
interaction region at all. Therefore we are not able to judge from 
our current simulations whether a detonation can form directly 
from a shock. Since typical detonation flames are of the order of 
a few meters in width and the best resolved regions in our sim- 
ulations are of the order of 10 5 m, only future simulations with 
very high adaptivity may be able to resolve their formation in 
global simulations. 

Secondly, it may be possible to form a detonation sponta- 
neously without a preceding shock. In this case, as first proposed 
by lZel'dovich et aT7| ( |197Q] >, nuclear burning in a preconditioned 



region (i.e. with a spatial temperature or composition gradient 
that causes an induction time gradient) creates a spontaneous ig- 
nition wave. If the preconditioning is favorable, and this ignition 
wave has a supersonic phase velocity, it forms a shock wave that 
is followed by nuclear burning. But only if the release of nu- 
clear energy behind the shock wave is sufficient to sustain the 
shock after it reaches the bottom of the gradient, it becomes a 
self-sustained detonation. 

Because the mechanism itself can not be resolved in our sim- 
ulations, we have to rely on studies that model the mechanism 
for different conditions in a microscopic region. Such studies 
provide an indication whether a detonation is expected to form 
for the physical conditions in our simulation. In a recent study, 
|Seitenzahl et al.| ((2009) carried out resolved one-dimensional 
simulations of the gradient mechanism for different initial con- 
ditions (i.e. varying temperature, density and steepness of the 
gradient). They show that for the most favorable conditions of 
geometry and gradient the density threshold to form a detona- 
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Fig. 3. Resolution test of the detonation conditions of the merger of two white dwarfs of 0.90 M and 0.89 M Q with an initial period 
of 25s. The left panel shows a simulation with 2 ■ 10 6 particles, the simulation shown in the right panel contains 10 7 particles. Both 
panels show a density slice after the system has evolved for 32s. The black stars mark all particles with temperatures above 2 • 10 9 K. 



tion in carbon-oxygen material is 10 6 g cirT 3 . At this density, a 
temperature above 2.8 • 10 9 K is required over a range of some 
meters. For higher densities, the temperatures required drop. In 
the following, we require a hot-spot with a density larger than 
2 ■ 10 6 g cirT 3 and a temperature above 2.5 • 10 9 K for reaching 
the conclusion that a detonation forms. 



4.2. Conditions reached in the merger simulation 

Figure [2] shows density versus temperature of all particles for 
several snapshots from 60 s to 90 s after the start of the simula- 
tion in a scatter plot. This is the time-range when the secondary 
white dwarf is disrupted, its material is violently accreted onto 
the primary white dwarf and a hot-spot at the interface emerges. 
Afterwards, the temperature drops again as the hot material ex- 
pands. Temperatures exceeding 2- 10 9 K are reached, but only for 
a small number of particles. These particles are the result of the 
additional energy release by thermonuclear combustion, because 
they ignite carbon which is subsequently burned. In simulations 
using identical initial conditions but no nuclear network they do 
not exist. Some of them are also in dense enough regions to ful- 
fill the criterion we require to form a detonation. 

As the overall number of these particles is quite small, their 
exact properties have to be taken with caution. Therefore we 
have to check whether they show a barely resolved physical ef- 
fect or are just a numerical artifact. This can be tested by redoing 
the simulation with significantly higher resolution. 

The binary is not completely merged at the time the first hot 
spots form. Therefore, if a detonation occurs, it stops any further 
evolution of the merger and the merged object never reaches the 
final state of a white dwarf surrounded by a hot envelope found 
in previous simulations (e.g. Guerrero et al. 2004 Yoon et al. 
2057||Loren-Aguilar et al.|2009) . 



4.3. Convergence study 

For a convergence study, we choose a slightly different binary 
with masses of 0.90 M and 0.89 M for the two WDs, respec- 
tively. As this binary has a mass ratio much closer to one, it is 
expected to merge more violently and therefore it should be eas- 
ier to find extreme conditions in it. 

Figure [3] shows the density of two simulations of such a 
merger with 2 ■ 10 6 and 10 7 particles, respectively. It demon- 
strates excellent agreement in the density structure between both 
simulations. In both simulations the particles with temperatures 
exceeding 2 • 10 9 K are found at the same locations. The number 
of these particles, however, increases from 11 to 127 from the 
lower resolved to the higher resolved simulation. As the number 
of particles is only increased by a factor of five, the total mass at 
temperatures above 2 • 10 9 K is approximately twice as large in 
the high resolution simulation. Therefore we can conclude that 
these hot particles are real, rather than a numerical artifact. The 
high resolution simulation indicates that we underestimate the 
conditions in these hot-spots due to a lack of resolution in our 
standard setup. 



4.4. Additional effects on the detonation conditions 

There are two additional effects that may influence the formation 
of a detonation, but are basically ignored in our simulations. 

As discussed above, we do not resolve shocks at the inter- 
face very well. This effect can only be overcome by choosing a 
scheme that allows for refinement according to arbitrary criteria. 

Second, according to stellar evolution calculations there is a 
small helium shell on top of a carbon-oxygen white dwarf. For 
an isolated white dwarf around 0.9 Mq this shell is expected to 
contain * 10~ 3 to 10~ 2 M G of helium jlben & Tutukov||l985t 
It is known from the latest studies by |Seitenzahl et aL] (2009 ) 
that already a small mass fraction of helium mixed into carbon 
fuel significantly lowers the density required for the formation 
of a detonation. As the expected helium layer is located exactly 
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Fig. 4. Temperature slice and density-temperature scatter diagram for four mergers with the same primary white dwarf, but different 
secondary white dwarfs, that have mass ratios of 0.99, 0.9, 0.84 and 0.78. The rows show a temperature slice (on the left) and the 
distribution of all particles in temperature-density space (on the right) approximately at the time when there are the most favorable 
conditions for a detonation. Red/blue colors indicate high/low temperatures. 

in the region where the merger is most violent and the hot-spots form, it may facilitate a detonation at even lower densities in our 

scenario. 
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Fig. 5. Final composition of the ejecta of the the explosion of a merger of two 0.89 M Q white dwarfs. At this time, 100s after 
the explosion, the ejecta expand homologously. The panels show from the top to the bottom the abundances of carbon, oxygen, 
intermediate mass elements and iron group elements. The left and right panels show the system viewed face-on and edge-on, 
respectively. 
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In summary, there are reasons to believe that detonations 
may form for even a wider range of conditions than the ones 
we have assumed here. 



6. Comparison with observations 

6. 1 . Composition of the ejecta 



5. The impact of the mass ratio on the merger 

The violent merger scenario can only be realized frequently 
enough to account for a significant fraction of all type la su- 
pernovae if it also works for moderate mass differences between 
the two white dwarfs. Thus, we have to check how changing the 
mass ratio for binary systems with a fixed primary white dwarf 
mass affects the conditions during the merger. Only if mergers 
with a reasonable range also lead to the formation of a detona- 
tion, the scenario will be more than an exotic possibility. To test 
this, we compare four different mergers. All have the same pri- 
mary mass of 0.9 M , but different secondary masses of 0.89 M , 
0.81 Mq, 0.76 M and 0.7 M G which give mass ratios of 0.99, 
0.9, 0.84 and 0.78. In our simulations, they are set up with initial 
periods of 25 s, 33 s, 36 s, and 40 s, respectively. 

Figure |4] shows densities and temperatures of all particles of 
these simulations at the time when the conditions are most fa- 
vorable for a detonation as well as temperature slices through 
the centers of the binaries. Obviously, there are considerable dif- 
ferences between these three systems. Firstly, with decreasing 
mass ratio, the merger becomes less violent. While the mergers 
with mass ratios of 0.99, 0.9 and 0.84 produce several hot par- 
ticles that ignite carbon and reach temperatures above 2 • 10 9 K, 
this is not the case for the merger with the smallest mass ra- 
tio. It is still possible that resolving the interaction region better 
will also show hotter particles. However, at the moment it seems 
more likely that below a certain mass ratio the merger is just not 
violent enough to ignite a detonation. This suggests a limiting 
mass ratio for the violent merger scenario of around 0.8. 

Another difference is the dynamical effect of the merger on 
the primary white dwarf. In a nearly equal mass merger the pri- 
mary star is heavily distorted, very similar to the equal mass 
merger described in Pakmor et al. ( 2010[ ). For a mass ratio of 
0.9, the primary white dwarf remains unaffected in the center, 
but its surface is distorted. For the smaller mass ratios, it stays 
completely intact and cool and is surrounded by the material of 
the disrupted less massive companion. 

The effect of the merger on the primary white dwarf can also 
be seen in the right panels of Figure|4] In the q = 0.99 merger the 
central density of the remaining white dwarf is lower than in the 
other cases. As for these mergers the amount of 56 Ni produced 
depends sensitively on the central density of the remaining white 
dwarf, this leads to smaller 56 Ni masses in mergers with higher 
mass ratios. 

While this does not change the scenario fundamentally, it 
breaks the relation between mass of the primary white dwarf 
and the final 56 Ni mass of the explosion for dim explosions. For 
explosions of mergers of more massive white dwarfs, however, 
where most of the 56 Ni is produced in nuclear statistical equilib- 
rium, this will be a minor effect only. 

As shown in Figure[4j with decreasing mass ratio the merger 
becomes less violent and the density in the hot spots drops. 
This may make the formation of a detonation more difficult for 
smaller mass ratios. 

For the smallest mass ratio no violent merger occurs and we 
see the onset of a different merger regime in which an accretion 
disk forms around primary white dwarf. 




Fig. 6. Mass distribution in velocity space resulting from the 
thermonuclear explosion following a violent merger of two 
white dwarfs of 0.9 M . The standard Chandrasekhar-mass SN 
la model W7 is overplotted for comparison. 




Fig. 7. Angle averaged composition in velocity space (straight 
lines) of the result of the thermonuclear explosion following a 
violent merger of two white dwarfs of 0.9 M . It is compared to 
the results of an abundance tomography study of SN 2005bl by 
|Hachinger et al.| ( [20 09 ) (dashed lines). Note that only a qualita- 
tive match can be expected, as the exact abundances depend on 
the density profiles, where Hachinger et al.| ( |2009| l used some- 
what different assumptions. 



Figure [5] shows the ejecta composition of the thermonuclear 
explosion of the merger of two 0.89 M as described in |Pakmo7 
|et al.| ( |2010| l, which is close to the merger with a mass ratio of 
0.99 described above. 

The final composition contains 0.03 M of carbon, 0.54 M 
of oxygen, 1.05 M of intermediate-mass elements and 0.1 M 
of iron group elements. As the densities at the time the explo- 
sion happens do not exceed 2 x 10 7 g crrT 3 , nuclear burning does 
not reach nuclear statistical equilibrium anywhere in the object. 
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Instead, iron group elements are only produced by incomplete 
silicon burning. 

As the nuclear burning in the detonation takes place at low 
densities, electron captures are not efficient. Therefore the ini- 
tial electron fraction Y e of the unburned material is conserved 
throughout the nuclear burning. In case of Y e = 0.5 the iron 
group elements consist of almost pure 56 Ni. A lower initial Y e 
is equivalent to more neutron-rich material in the pre-explosion 
composition (i.e. more 22 Ne resulting from the He burning 
phase). It leads to the production of some amount of stable iron 
replacing some of the 56 Ni. In total, however, the fraction of sta- 
ble iron group elements produced is very small. 



As shown in Figure [5 
center only, surrounded " 



iron group elements are located at the 
->y intermediate-mass elements which 
are mixed with oxygen. Most carbon is found in the outermost 
parts, but there is some mixing with intermediate-mass elements. 
It is important to note that both, carbon and oxygen, are at dif- 
ferent locations than the iron group elements, i.e. there is no 
mixing of unburned material with iron group elements. There 
is also a clear difference between the distribution in the x-y- 
plane, which is the plane of rotation, and the distribution along 
the z-axis. Along the z-axis, the iron group elements are spread 
out much more, which means that in this direction there are 
iron group elements at higher velocities than within the x-y- 
plane. The opposite trend is observed for oxygen and interme- 
diate mass elements, which reach further out in the plane of 
rotation than perpendicular to it. This is a result of the den- 
sity structure of the merged object through which the detona- 
tion propagates. Consequently, lightcurves and spectra of this 
explosion are expected to show considerable viewing angle de- 
pendence (for viewing angle dependent lightcurves see (Pakmor 
|et al.|2010) . 



6.2. Comparison with other models and observations 



Compared to the standard W7 explosion model of Nomoto et al. 
(1984) our explosion models contains about 25% more mass, 
but the kinetic energy of the ejecta is about the same. Therefore 
the velocities are lower on average. Figure [6] shows the mass 
distribution in velocity space in comparison to the W7 model. 
The merger has more mass at velocities below 12, 000km s -1 
and significantly less mass at higher velocities. With a 56 Ni 
yield of only around 0.1 M Q it is in the range of subluminous 
SNe la. We argue that it falls into the class of 1991bg-like su- 
pernovae de scribed by |Leibundg ut et al.| ( |1993| l; |Taub enberger 
et~alj ( [2008"l e.g.). For one 1991bg-like super nova, SN 2005M, 
a detailed spectral analysis was presented by Hachinger et al. 
(2009). In this study the authors find that decent models for the 
observed spectra require significantly less mass at higher veloc- 
ities than what is inferred from the W7 mass profile. This is in 
agreement with the results of our simulations. Figure [6] shows a 
comparison between our final abundance pattern and the results 
of the abundance tomography for SN 2005bl. Despite differences 
in details, the main features are in qualitative agreement. There 
are basically three layers that differ in composition: The out- 
ermost layers are dominated by oxygen, with contributions of 
unburned carbon and some amount of intermediate mass ele- 
ments. Below about 11,000km s _I intermediate mass elements 
dominate down to very low velocities. Carbon is found down to 
* 9,000 km s 1 . Iron group elements are present up to veloci- 
ties of m 7,500km s _1 . In our model they dominate only in the 
very innermost part of the ejecta that is not probed by the photo- 
spheric spectra used in |Hachinger et al.] ( |2009| ). 



The composition pattern of the ejecta of SN 2005bl inferred 
from abundance tomography is suggestive of a detonation at low 
densities as provided by our model. In particular, the fact that 
iron group elements co-exist with intermediate-mass elements 
down to very low velocities but - at least as far as the obser- 
vations can probe it - never dominate would be hard to explain 
with mixing alone. It strongly points to an extended range of 
incomplete silicon burning. 

Note that the abundances from our simulations shown in 
Fig. [6] are angle averaged. Taking into account the asymmetry 
of the explosion as shown in Fig [5] there is obviously some 
angle dependence in the composition in velocity space. This 
asymmetry will also affect the polarization of the observed spec- 
tra. Even though no spectrapolarimetry data is available from 
SN 1991bg or SN 2005bl there exist polarization spectra from 
SN 1999by, which falls into the same class of objects. |Howell| 
|et al.| ( |2001| ) find that this supernova shows continuum polar- 
ization of « 0.3% - 0.8%. This marks an important difference 
to normal SNe la, which do not show continuum polarization 
(Wang & Wheeler 2008 ). Our model provides a natural explana- 
tion for the comparatively high degree of polarization in spectra 
of 1991bg-like objects, as the explosion is intrinsically asym- 
metric. 

As this model assumes a progenitor system that is most likely 
different from that of normal SNe la, it may also be able to ex- 
plain the observed trend of 1991bg-like SNe to favor early-type 
galaxies and therefore old stellar populations (e.g. Howell 2001[ ). 
A long inspiral time of the original WD-WD binary leads to a 
long delay time between the formation of the progenitor system 
and the supernova explosion in our model. 

Overall, the violent merger model presented here is not only 
the first realistic hydrodynamical model of a thermonuclear ex- 
plosion of a WD that is able to produce 56 Ni masses small 
enough to explain subluminous SN la, but it also shows remark- 
able agreement with observable features of 1991bg-like objects 
in a more detailed comparison. A more quantitative compari- 
son based on synthetic lightcurves and spectra can be found in 
|Pakmor et d.] ( |2010) . 



7. Conclusion 

In this paper we investigate the effect of variations in the mass 
ratio of merging white dwarf binaries on the violent merger sce- 
nario proposed by Pakm or et al.| ( [20 10 ). We first concentrate on a 
prototypical realization of this scenario: the case of a double de- 
generate binary system with a total mass of 1.61 M Q and a mass 
ratio of 0.9. We follow the evolution of the binary system and 
subsequent merging event in smoothed particle hydrodynamical 
simulations. 

We discuss in detail the conditions that arise in the violent 
merging event and find that it is reasonable to assume that they 
lead to the formation of a detonation. We also show that while 
we are only barely resolving the most important part of the sim- 
ulation (i.e. the interface between the two white dwarfs, where 
the detonation forms), a very high resolution run comprising 10 7 
particles turns out to be even more in favor of the formation of a 
detonation. 

To estimate the influence on the mass ratio of the merging 
WDs on the detonation ignition conditions, we compare mergers 
with a primary white dwarf mass of 0.9 M , but different mass 
ratios ranging from 0.76 to unity. We find that merging systems 
with a mass ratio above ~ 0.8 detonate, but those with a mass 
ratio below 0.8 fail to do so. 
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Finally, we compare the final abundance structure in the 
ejecta of a thermonuclear explosion after the violent merger of 
two white dwarfs with ~ 0.9 M Q with observed properties of 
faint SN 1991bg-like supernovae. The ejecta have higher mass, 
but smaller velocities compared to the standard Chandrasekhar- 
mass W7 model. The abundance distribution of the model pro- 
vides a good qualitative match with the composition found from 



the abundance tomography of SN 2005bl by |Hachinger et al. 



(2009). Moreover, the model could provide a natural explana- 
tion for the polarization found in SN 1999by and connect the 
preference of SN 1991bg-like objects to occur in old stellar pop- 
ulations with their progenitors. 

Of course, future studies are needed to explore the parameter 
space of different white dwarf masses and mass ratios in the vi- 
olent merger scenario. Only this will allow us to understand the 
complete range of possible outcomes and the contribution of this 
channel to the observed sample of SNe la. 
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